Correlations and Pair Formation in a Repulsively Interacting Fermi Gas 
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A degenerate Fermi gas is rapidly quenched into the regime of strong effective repulsion near a 
Feshbach resonance. The spin fluctuations are monitored using speckle imaging and, contrary to 
several theoretical predictions, the samples remain in the paramagnetic phase for arbitrarily large 
scattering length. Over a wide range of interaction strengths a rapid decay into bound pairs is 
observed over times on the order of f %/Ef, preventing the study of equilibrium phases of strongly 
repulsive fermions. Our work suggests that a Fermi gas with strong short-range repulsive interactions 
does not undergo a ferromagnetic phase transition. 
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Many-body systems can often be modeled using con- 
tact interactions, greatly simplifying the analysis while 
maintaining the essence of the phenomenon to be studied. 
Such models are almost exactly realized with ultracold 
gases due to the large ratio of the de Broglie wavelength 
to the range of the interatomic forces pQ. For itinerant 
fermions with strong short-range repulsion, textbook cal- 
culations predict a ferromagnetic phase transition - the 
so-called Stoner instability [2 . 

Here we investigate this system using an ultracold gas 
of fermionic lithium atoms, and observe that the ferro- 
magnetic phase transition does not occur. A previous 
experimental study [3 employing a different apparatus 
found indirect evidence for a ferromagnetic phase, but 
did not observe the expected domain structure, possibly 
due to the lack of imaging resolution. Here we address 
this shortcoming by analyzing density and spin density 
fluctuations via speckle imaging 4|. When spin domains 
of m atoms form, the spin density variance will increase 
by a factor of m [S], even if individual domains are not 
resolved. One main result of this paper is the absence 
of such a significant increase which seems to exclude the 
possibility of a ferromagnetic state in the studied system. 

The Stoner model assumes a two-component Fermi gas 
with a repulsive short-range interaction described by a 
single parameter, the scattering length. The predicted 
phase transition to a ferromagnetic states requires large 
repulsive scattering lengths on the order of the inter- 
atomic spacing. They can be realized only by short- 
range attractive potentials with a loosely bound state 
with binding energy ft 2 /(ma 2 ), with m being the atomic 
mass and a being the scattering length [BJ. However, the 
repulsive gas is then by necessity only metastable with 
respect to decay into the bound state. Many theoreti- 
cal studies of a Fermi gas with strong short-range repul- 
sive interactions assume that the metastable state is suf- 
ficiently long-lived [7HT5] . In recent Monte-Carlo simula- 
tions, the paired state is projected out in the time evolu- 
tion of the system 19, 20J . Theoretical studies concluded 
that the pairing instability is somewhat faster than the 
ferromagnetic instability [2"Tl |2"2] . The second major re- 



sult of this paper is to show that pair formation occurs 
indeed on a very short time scale. The measured time 
constant of lOh/Ep (where Ep is the Fermi energy) indi- 
cates that the metastable repulsive state will never reach 
equilibrium and that, even in a metastable sense, a Fermi 
gas with strong short-range repulsive interactions does 
not exist. The fast pair formation could not be observed 
previously due to limited time resolution |3]. Instead, a 
much slower second phase in the conversion of atoms to 
pairs was observed leading to the wrong conclusion that 
the unpaired atoms have a much longer lifetime. 
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FIG. 1: Diagram showing energy levels and timing of the 
experiment. The upper (repulsive) and lower (attractive) 
branch energies, near a Feshbach resonance, are connected by 
3-body collisions. In our experiment, we quickly jump from 
a weakly interacting Fermi gas (A) to a strongly interacting 
one (B) with a rapid magnetic held change. The evolution of 
correlations and domains and the molecule formation (popu- 
lation of the lower branch) are studied as a function of hold 
time t. Adapted from [23] . 

The experiments were carried out with typically 4.2 x 
10 5 6 Li atoms in each of the two lower spin states |1) 
and 1 2) confined in an optical dipole trap with radial 
and axial trap frequencies tu r — 2ir x 100(l)s and 
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uj z = 2n x 9.06(25)s _1 . The sample was evaporatively 
cooled at a magnetic bias field B = 320G, identical 
to the procedure described in [24]. Then the magnetic 
field was slowly ramped to 730 G (k F a = 0.35) in 500 
ms. The fraction of atoms being converted to molecules 
during the ramp was measured (see below for method) 
to be below 5 %. The temperature of the cloud was 
typically 0.23(3) Tp at 527 G with a Fermi energy of 
Ep = kgTp = h x 6.1kHz. After rapidly switching the 
magnetic field from 730 G to the final value in less than 
350 /is, spin fluctuations were measured by speckle imag- 
ing. Optionally an appropriate RF pulse was applied di- 
rectly before imaging to rotate the spin orientation along 
the measurement axis. Due to the use of 20 cm diameter 
coils outside the vacuum chamber, the inductance of the 
magnet coils was 330 /iH and the fast switching was ac- 
complished by rapidly discharging capacitors charged to 
500V. 

Experimentally, spin fluctuations are measured using 
the technique of speckle imaging described in Ref. [4j. 
For an appropriate choice of detuning, an incident laser 
beam experiences a shift of the refractive index propor- 
tional to the difference between the local populations of 
the two spin states N% and N%. Spin fluctuations create 
spatial fluctuations in the local refractive index and im- 
print a phase pattern into the incoming light, which is 
then converted into an amplitude pattern during prop- 
agation. The resulting spatial fluctuations in the probe 
laser intensity are used to determine the spin fluctuations 
in the sample. 

In Ref. [4] we prepared samples on the lower branch 
of the Feshbach resonance, where positive values of kpa 
correspond to a gas of weakly bound molecules. At 
kpa = 1.2, we observed a sixfold suppression of spin 
fluctuations and a fourfold enhancement of density fluc- 
tuations. Typical fluctuations in the speckle images of 
a non-interacting Fermi gas at T — 0.23Tp amount to 5 
% of the average optical signal per pixel, corresponding 
to about 50 % of Poissonian fluctuations. Those fluctua- 
tions are modified by factors between 0.2 and 1.6 due to 
pairing and interactions. 

In this study, on the upper branch of the Fesh- 
bach resonance, the situation is reversed. For unbound 
atoms, as the interaction strength increases, the two 
spin components should develop stronger and stronger 
anticorrelations and enhanced spin fluctuations. Previ- 
ous experimental work [3] and several theoretical studies 
[TU1 [TT1 H3HE3 US] predicted a phase transition to a fer- 
romagnetic state where the magnetic susceptibility and 
therefore the spin fluctuations diverge. Recent Monte 
Carlo simulations [H] predict such a divergence around 
kpa = 0.83. We therefore expected an increase of spin 
fluctuations by one or several orders of magnitude, re- 
lated to the size of magnetic domains. 

Figure [2] shows the observed spin fluctuations enhance- 
ment compared to the non-interacting cloud at 527 G. 
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FIG. 2: Spin fluctuations (a) after 350 /j,s as a function of 
magnetic field and (b) on resonance as a function of hold time 
scaled to the value measured at 527G. Even at strong repul- 
sive interactions, the measured spin fluctuations are barely 
enhanced, indicating only short-range correlations and no do- 
main formation. The spin fluctuations were determined for 
square bins of 2.6 fim, each containing on average 1000 atoms 
per spin state. 



The variance enhancement factor reaches its maximum 
value of 1.6 immediately after the quench, decreasing 
during the 2 ms afterward. The absence of a dramatic 
increase shoes that no domains form and that the sample 
remains in the paramagnetic phase throughout. Similar 
observations were made for a wide range of interaction 
strengths and wait times. Note that first-order perturba- 
tion theory [35] predicts an increase of the susceptibility 
by a factor of 1.5 at kpa = 0.5 and by a factor of 2 
at kFa — 0.8 (i.e. no dramatic increase for kpa < 1). 
Therefore, our data shows no evidence for the Fermi gas 
approaching the Stoner instability. 

Before we can fully interpret these findings, we have 
to take into account the decay of the atomic sample on 
the upper branch of the Feshbach resonance into bound 
pairs. We characterize the pair formation by compar- 
ing the total number of atoms and molecules N a +2N mo \ 
(determined by taking an absorption image after bal- 
listic expansion at high magnetic field where molecules 
and atoms have the same absorption resonance) to the 
number of free atoms (determined by rapidly sweeping 
the magnetic field to 5G before releasing the atoms and 
imaging the cloud, converting pairs into deeply bound 
molecules that are completely shifted out of resonance) 

The time evolution of the molecule production (Fig. 
[3]) shows two regimes of distinct behavior. For times 
less than 1 ms, we observe a considerable number of 
atoms converted into molecules, while the total number 
N a +2N rno \ remains constant. The initial drop in atom 
number becomes larger as we increase the final magnetic 
field, and saturates at around 50 % near the Feshbach 
resonance. 
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We attribute this fast initial decay in atom number to 
recombination [2JJ [28] into the weakly bound molecular 
state. We obtain an atom loss rate N a /N a = 250 s _1 
at 790G in the first lms after the magnetic field switch. 
Assuming a three-body process we estimate the rate co- 
efficient L 3 at this field to be 3.9 x 10~ 22 cm 6 s _1 , though 
the interaction is already sufficiently strong for many- 
body effects to be significant. For stronger interactions, 
about 30% of atom loss occurs already during the rele- 
vant 100 (is of ramping through the strongly interacting 
region, indicating a lower bound of around 3 x 10 3 s _1 
for the loss rate which is 13% of the inverse Fermi time 
Ep/Ti, calculated with a cloud averaged Fermi energy. 

After the first millisecond, the molecule formation rate 
slows down, by an order of magnitude at a magnetic field 
of 790G (and even more dramatically at higher fields) 
when it reaches about 50 %. It seems likely that the 
molecule fraction has reached a quasi-equilibrium value 
at the local temperature, which is larger than the ini- 
tial temperature due to local heating accompanying the 
molecule formation. Ref. [29] presents a simple model for 
the equilibrium between atoms and molecules (ignoring 
strong interactions). For phase space densities around 
unity and close to resonance, the predicted molecule frac- 
tion is 0.5, in good agreement with our observations [50] , 

For longer time scales (hundred milliseconds) we ob- 
serve a steady increase of the molecule fraction to 90 % 
for the longest hold time. This occurs due to continuous 
evaporation which cools down the system and shifts the 
atom-molecule equilibrium towards high molecule frac- 
tions. During the same time scale, a slow loss in both 
atom number and total number is observed caused by 
inelastic collisions (vibrational relaxation of molecules) 
and evaporation loss. 

Is the rapid conversion into molecules necessarily faster 
than the evolution of ferromagnetic domains? Our an- 
swer is tentatively yes. First, for strong interactions with 
kpa around 1, one expects both instabilities (pair forma- 
tion and Stoner instability) to have rates which scale with 
the Fermi energy Ep and therefore with n 2 / 3 . Therefore, 
one cannot change the competition between the instabil- 
ities by working at higher or lower densities. According 
to Ref. [3T] the fastest unstable modes for domain for- 
mation have a wavevector ga and grow at a rate of 
up to Ep/Ah when the cloud is quenched sufficiently far 
beyond the critical interaction strength. Unstable modes 
with such wavevectors will develop "domains" of half a 
wavelength or size £ = n/q = 2tt /kp containing 5 atoms 
per spin state in a volume £ 3 . This rate is comparable 
to the observed conversion rates into pairs of O.l'SEp. 
Therefore, at best, "domains" of a few particles could 
form, but before they can grow further and prevent the 
formation of pairs (in a fully polarized state), rapid pair 
formation takes over and populates the lower branch of 
the Feshbach resonance. Based on our observations and 
these arguments, it seems that it is not possible to realize 
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FIG. 3: Characterization of molecule formation at short and 
long hold times, and at different values of the interaction 
strength. The closed symbols, circles (black) at 790G with 
kpa = 1.14, squares (blue) at 810G with kpa — 2.27 and dia- 
monds (red) at 818G with kpa — 3.5 represent the normalized 
number of free atoms, the open symbols the total number of 
atoms including those bound in Feshbach molecules (open cir- 
cles at 790G with kpa = 1.14). The crosses (green) show the 
molecule fraction. The characteristic time scale is set by the 
Fermi time h/Ep = 43fis, calculated with a cloud averaged 
Fermi energy. 



ferromagnetism with strong short range interaction, and 
therefore the basic Stoner model cannot be realized in 
nature. 

One possibility to suppress pair formation is provided 
by narrow Feshbach resonances. Here the pairs have 
dominantly closed channel character and therefore a 
much smaller overlap matrix element with the free atoms. 
However, narrow Feshbach resonances are characterized 
by a long effective range and do not realize the Stoner 
model which assumes short-range interactions. Other 
interesting topics for future research on ferromagnetism 
and pair formation include the effects of dimensionality 
[3T1 I3"2"] . spin imbalance [33] 131], mass imbalance [55] . 
lattice and band structure [35] [37] . 

We now discuss whether ferromagnetism is possible af- 
ter atoms and molecules have rapidly established local 
equilibrium. In other words, starting at T = 0, one could 
heat up the fully paired and superfluid system and cre- 
ate a gas of atomic quasiparticles which are similar to 
free atoms with repulsive interactions. Density and tem- 
perature of the atoms are now coupled. It is likely that 
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such a state is realized in our experiments after a few 
ms following the quench, until evaporative cooling con- 
verts the system into a molecular condensate over sa 100 
ms. The possibility that such a quasiparticle gas could 
become ferromagnetic has not been discussed in the liter- 
ature. Our experiments do not reveal any major increase 
in spin fluctuations which seems to exclude a ferromag- 
netic state. In the simplest picture, we could regard the 
atomic quasiparticles as free atoms, and then apply the 
Stoner model to them. Ferromagnetic domain formation 
is analogous to phase separation between the two spin 
components [3J. Since dimers interact equally with the 
two spin components, one might expect that even a no- 
ticeable dimer fraction should not suppress the tendency 
of the atomic gas to form domains. Therefore, in a simple 
model, one may neglect dimer-atom interactions. 

If the Stoner model applies to this quasiparticle gas, 
the next question is whether the temperature is low 
enough for the ferromagnetic phase transition. Avail- 
able theoretical treatments do not predict an exact max- 
imum transition temperature to the ferromagnetic state 
and obtain an unphysical divergence for large scatter- 
ing lengths. Since the only energy scale is the Fermi 
temperature, one would expect a transition temperature 
which is a fraction of the Fermi temperature |38] , higher 
or around the temperature scale probed in our experi- 
ments. However, even above the transition temperature, 
the susceptibility is enhanced. A simple Weiss mean field 
or Stoner model leads to the generic form of the suscep- 
tibility X (T) = Xo(T)/(l ~ wXo{T)), where Xo (T) is the 
Pauli susceptibility of the non-interacting gas and w the 
interaction parameter. This formula predicts a two-fold 
increase in the susceptibility even 50 % above the transi- 
tion temperature, which is well within the sensitivity of 
our measurements. 

Therefore, our experiment can rule out ferromagnetism 
for temperatures even slightly lower than the experimen- 
tal temperatures. Temperatures are very difficult to mea- 
sure in a transient way for a dynamic system which may 
not be in full equilibrium. For example, cloud ther- 
mometry requires full equilibration and lifetimes much 
longer than the longest trap period. We attempted to 
measure the temperature after the hold time near the 
Feshbach resonance by quickly switching the magnetic 
field to weak interactions at 527 G and then performing 
noise thermometry using speckle imaging [?]. We mea- 
sure column-integrated fluctuations that are 0.61(8) of 
the Poisson value which implies an effective temperature 
well below Tp, around 0.33(7) Tp, not much higher than 
our initial temperature of 0.23 Tp. Although the cloud 
is not in full equilibrium, an effective local temperature 
can still be obtained from noise thermometry. 

Alternatively, we can estimate the temperature in- 
crease from the heat released by pair formation. A sim- 
ple model [39j accounting for the relevant energy con- 
tributions predicts for kpa — 1 that molecule fractions 



of higher than 20% result in a final temperature above 
0.4Tp, an estimate which is higher than the measurement 
reported above. One may hope that closer to resonance 
many-body effects lower the released energy, however as 
we show in the supplement (Fig. 1 of [3S]) this is not 
necessarily the case due to the repulsive interaction en- 
ergy. 

Our experiment has not shown any evidence for a pos- 
sible ferromagnetic phase in an atomic gas in "chemical" 
equilibrium with dimers. This implies one of the follow- 
ing possibilities, (i) This gas can be described by a simple 
Hamiltoninan with strong short range repulsion. How- 
ever, this Hamiltonian does not lead to ferromagnetism. 
This would be in conflict with the results of recent quan- 
tum Monte-Carlo simulations [TH1 [20] and second order 
perturbation theory |llj . and in agreement with conclu- 
sions based on Tan relations [3D], (ii) The temperature 
of the gas was too high to observe ferromagnetism. This 
would then imply a critical temperature around or below 
0.2T/Tp, lower than generally assumed, (iii) The quasi- 
particles cannot be described by the simple model of an 
atomic gas with short-range repulsive interactions due to 
their interactions with the paired fraction. 

A previous experiment [3] reported evidence for fer- 
romagnetism by presenting non-monotonic behavior of 
atom loss rate, kinetic energy and cloud size when ap- 
proaching the Feshbach resonance, in agreement with 
predictions based on the Stoner model. Our measure- 
ments confirm that the properties of the gas strongly 
change near kpa = 1. Similar to [3], we observe features 
in kinetic and release energy measurements near the reso- 
nance (see Supplemental Information 39J). However, the 
behavior is more complex than that captured by simple 
models. The atomic fraction decays non-exponentially 
(see Fig. [3]), and therefore an extracted decay time will 
depend on the details of the measurement such as time 
resolution. Ref. [3 found a maximum of the loss rate of 
200 s _1 for a Fermi energy of 28 kHz. Our lower bound 
of the decay rate of 3 x 10 3 s _1 is 15 times faster at a five 
times smaller Fermi energy. Our more detailed study 
rules out that Ref. [3] has observed ferromagnetic be- 
havior. 

Our conclusion is that an ultracold gas with strong 
short range repulsive interactions near a Feshbach res- 
onance remains in the paramagnetic phase. The fast 
formation of molecules and the accompanying heating 
makes it impossible to study such a gas in equilibrium, 
confirming predictions of a rapid conversion of the atomic 
gas to pairs [HIST] . Therefore, it appears that the widely 
used Stoner model cannot be realized in Nature in its sim- 
plest form since the neglected competition with pairing 
is crucial. 
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